On Dark Matter Annihilation in the Local Group. 



cn 
o 
o 

> 

o 



(N 
> 

o . 

(N ■ 
l> ■ 

o ; 
m . 
o 

■ 

Oh; 

6 ■ 



Lidia Pieri^'^'0 and Enzo Branchini^ 

'Dept. of Physics, Universita di Roma Tre, Rome, Raly 
^INFN Roma Tre, Rome, Raly 

Under the hypothesis of a Dark Matter composed by supersymmetric particles like neutralinos, 
we investigate the possibility that their annihilation in the halos of nearby galaxies could produce 
detectable fluxes of 7-photons. Expected fluxes depend on several, poorly known quantities such as 
the density proflles of Dark Matter halos, the existence and prominence of central density cusps and 
the presence of a population of sub-halos. We find that, for all reasonable choices of Dark Matter 
halo models, the intensity of the 7-ray flux from some of the nearest extragalactic objects, like M31, 
is comparable or higher than the difi^use Galactic foreground. We show that next generation ground- 
based experiments could have the sensitivity to reveal such fluxes which could help us unveiling the 
nature of Dark Matter particles. 



I. INTRODUCTION 



The nature of the Dark Matter (DM) is a fundamen- 
tal missing piece of the dark puzzle of the universe, and 
represents one of the hottest challenges facing particle 
physics and cosmology today. The latest observational 
data prefer a flat universe with a 23% mass fraction 
of DM, seeded with Gaussian, scale dependent adiabatic 
perturbations. The bulk of DM is believed to be "cold" , 
i.e. composed by weakly interacting particles that were 
non-relativistic at the epoch of decoupling. The most 
popular candidate for Cold Dark Matter (CDM) is the 
lightest supersymmetric particle (LSP) which, in most 
supersymmetry breaking scenarios, is the neutralino x- 
Neutralinos are spin-i Major ana fermions, linear com- 



bination of the neutral gauge bosons and neutral Higgs 
doublet spartners, x° = aB + bW^ + cH^ + dH^. In the 
popular SUGRA or SUGRA-like models, where gaugino- 
universality is required, its mass is constrained by accel- 
erator searches and by theoretical considerations of ther- 
mal freeze-out to lie in the range 50 GeV < < 10 TeV 
0, 1^. Neutralinos decouple at a temperature that is 
roughly hence they behave like CDM. If i?-parity 

is conserved, neutralinos are stable and can change their 
cosmological abundance only through annihilation, which 
implies that the resulting production of detectable mat- 
ter, antimatter, neutrinos and photons is enhanced in 
high density regions like the center of virialized DM halos. 
At low energies, direct detection terrestrial experiments 
are being performed to measure the energy deposited by 
elastic neutralino- nucleon scattering 13, IBl. Besides di- 
rect searches, indirect detection of DM through its an- 
nihilation products can be implemented 0, Q, which 
are necessary for energies greater than 250 GeV. Most 
indirect searches rely on detection of antimatter | a Is l 
and 7-rays produced in the Milky Way halo 0, llll Haj 
or in Galactic sub-structures [131 [iJI ; the possibility of 
looking at extragalactic sources has also been suggested 
[li lli lid lli llM llSl 12^- Unfortunately, our partial 
knowledge of those astrophysical processes that govern 



the assembly of DM concentrations did not allow to make 
firm predictions on DM detectability. In particular, the 
central structure of the DM halos and the prominence of 
small scale structures within the smooth DM halo are far 
from being well determined. The purpose of this work 
is to explore the possibility of extending DM searches 
to extragalactic sources in our Local Group of galaxies 
(LG) and to make robust predictions about their actual 
detectability with ground-based experiments once model 
uncertainties are taken into account. Supersymmetric 
LSP are not the only particles which can produce anni- 
hilation signals (e.g. |21|), yet in this work we focus on 
the contribution of neutralinos only. 

The outline of the paper is as follows: in section II 
we discuss the theoretical framework and present the 
predictions from the smooth halo models. In section 
III and IV we show how predictions change when the 
presence of Supermassive Black Holes (SMBH) and 
that of a population sub-galactic halos are taken into 
account. Section V deals with the actual detectability 
of LG sources and Galactic sub-halos by ground-based 
detectors. Our conclusions are presented in Section VI. 



II. THEORETICAL FRAMEWORK 

In a phenomenological approach which considers the LSP 
as a DM candidate, the expected photon flux from neu- 
tralino annihilation is given by: 
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In this formula, the dependence of the flux on particle 
physics inputs and on cosmological inputs are factorized. 
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Particle physics. The first sum in brackets represents 
a 7-line, i.e. the 1-loop mediated process XX ~^ 77 or 
YY — ^ Zj. Since neutraUnos move at Galactic speed 
[22], their annihilation occurs at rest and the outgoing 
photons carry an energy equal (77 final state) or very 
close {Z'y final state) to the neutralino mass. For the 
corresponding 7-lines we have N^j = 2 and Nz-y = 1 
monochromatic final state photons. Branching ratios 
lO^'^ imply that the photon flux is dom- 
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inated by the continuum emission rather than by the 7- 
lines. The continuum emission is given by the second 
sum of Eq. |^ running over all the tree-level final states. 
At the tree- level neutralinos annihilate into fermions, 
gauge bosons, Higgs particles and gluons. Decay and/or 
hadronization in tt'^ give a continuum spectrum of 7- 
photons emerging from the tt^ decay. Depending on the 
annihilation channel, the continuum differential energy 
spectrum can be parametrized as dNj/dx 
where x — E-y/m 



X, a ~ 0(1) and b - 0(10). These 
spectra have been calculated in using a PYTHIA 
Monte Carlo simulation. Branching ratio bp calculations 
for each process involved in the sum and the evaluation 
of the total annihilation cross section (aav) depend on 
the assumed supersymmetric model. Here we consider 
as a unique annihilation channel the one into W bosons, 
XX ~^ W'^W' , with parameters a = 0.73 and b = 7.76 in 
the above mentioned dN^/dx formula. We also use the 
constraint to the value of the thermally-averaged annihi- 
lation cross-section (ctow) given by the neutralino number 
density at the freeze-out epoch jj]: 
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where h is the Hubble constant, H, in units of 
100 km Mpc~^ and fi^ is the average density of 
neutralinos in units of critical density, pc = /SttG. 
Unless otherwise specified, throughout the paper we 
consider neutralinos with a mass of 1 TeV and cross 
section (ctqw) — 2 • 10^^^ cm^^ s^^ . 

Cosmology. The last term in Eq. (Q) describes the ge- 
ometry of the problem and assumes that the DM is con- 
centrated in a single, spherical DM halo of radius r^ax 
and density profile Px(^) located at distance D from the 
observer. The integral can be easily generalized to any 
DM distribution once that the DM density along the line 
of sight, Px{'>')i is specified. 

An accurate evaluation of this integral should account 
for the fact that DM halos are triaxial rather than spheri- 
cal objects 122] . Deviations from sphericity do not change 
appreciably the 7-ray flux from LG objects within the 
typical acceptance angle of the detectors (0.1° — 1°) and 
would only enhance the Galactic flux by a modest 15% 
[i^ l and thus will not be considered here. The shapes of 
the halo density profiles are poorly constrained by obser- 
vations, especially near the center. Indeed, recent mea- 



surements of the dwarf galaxy rotation curves, as well 
as of the intracluster medium and of the gravitational 
lensing, could not discriminate between dark halos with 
constant density core and cusps [l^. Similarly, ana- 
lytical arguments did not yet lead to a definitive, unique 
prediction [illiTlEit. 

The most stringent constraints on the shape of DM 
halos are currently obtained from N-body simulations. 
Although the resolution in numerical experiments is still 
an issue |29l |. advances in numerical techniques have al- 
lowed to discover that DM halos have universal density 
shapes p^if) Ps{f/fs)~"- While singular isothermal 
halos characterized by a steep cusp are ruled out, cur- 
rent numerical experiments do not allow to discriminate 
among the Moore profile (3^] characterized by a = 1.5 
and the shallower NEW profile [sJl with a = 1 (or slightly 
shallower ones .24]). In this paper we assume that the 
actual shape of the DM density profiles is bracketed by 
the Moore and the NEW models and thus we can obtain 
a fair estimate of current model uncertainties by consid- 
ering both of them in our analysis. 

Both models have central cusps where the neutralino an- 
nihilation rate is greatly enhanced. The scale radius, 
rg, and the scale density, ps, can be fixed by observa- 
tions (the virial mass of the halo or its rotation velocity) 
and theoretical considerations that allow to determine 
the concentration parameter c = Tmr/fs (where the virial 
radius, r^ir, is defined as the radius within which the halo 
average density is 200pc)- 

Since the largest 7-ray fluxes are expected from nearby 
objects, we have considered only the 44 nearest LG galax- 
ies and their parent halos. Table HI lists the masses, posi- 
tions (from -32]) and virial radii of the galaxies which are 
more relevant for our analysis, along with their scale radii 
and scale densities computed assuming a Moore profile. 
Their concentration parameters, cnfw and C Mnnr p = 
0.64 cnfWi have been computed according to [2^1 as- 
suming CDM power spectrum with a shape parameter 
r = 0.2, normalized to ug = 0.9. In our analysis we have 
also included the giant elliptical galaxy M87 at the cen- 
ter of the Virgo cluster whose predicted annihilation flux 
was already computed by [l9j |. The mass, distance and 
the virial radius of M87, listed in Tabled are taken from 
0. 

Since the intensity of the 7-ray flux is proportional to 
p^, great attention must be paid in modeling the central 
region of the halo. The fact that the Moore model pre- 
dicts a divergent flux from the halo center implies that 
there must exist a minimum radius, r„ 



within which 
- {{(Jav) n^irmin)y^ cquals 



the self-annihilation rate ti 

the dynamical time tdy,i ~ (Gpv)~i , where p-^ is the 
mean halo density within rmin |35| . We have used this 
prescription to determine rmin for all our 44 halos and 
model their DM distribution with a NEW or a Moore 
profile with a small constant core Pxi'i'min)- The result- 
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FIG. 1: Aitoff projection of the 44 nearest LG galaxies in 
Galactic coordinates. The size of each symbol is scaled to the 
7-ray flux emitted by a host DM halo with a Moore profile 
within a viewing angle of 1° from the halo center. 



ing profiles are 
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for the NFW case, and 
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Fig. 121 shows the smooth 7-ray Galactic foreground 
for a Moore (continuous line) and a NFW (dotted line) 
profile as a function of ip, the angle from the Galactic 
Center (GC). Filled triangles represent the 7-fluxes from 
those LG galaxies which are brighter than the Galac- 
tic foreground and from M87, assuming a Moore profile. 
The fluxes from the Small and Large Magellanic Clouds 
(SMC, LMC), M31 and, to a lesser extent, M87 are above 
the Galactic level. In some cases, like M33 and Sagittar- 
ius, extragalactic and Galactic contributions are compa- 
rable. Similar considerations are valid for the case of a 
NFW profile (filled dots), despite of the significant re- 
duction of the extragalactic fluxes that cause M33 and 
Sagittarius to be dimmer than the Galactic foreground. 
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49 
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6.8 
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M31 


2.0 • 10^^ 


770 


258 


47.3 


8.6 ■ 10^ 


M33 


4.0 • 10^° 


840 


70 


10.6 


1.6 ■ 10^ 


M87 


4.2 • 10" 


1.5 ■ 10^ 


1.55 ■ 10^ 


461 


2.6 ■ 10^ 



TABLE 1: Input parameters for the 7-ray flux prediction are 
shown for our Galaxy, M87 and for the LG galaxies brighter 
than the Galactic foreground. Masses and distances are taken 
from [3^. Virial radii were calculated following the prescrip- 
tion given in the text. Scale radii and scale densities were 
computed assuming a Moore profile. 



III. BLACK HOLES, CUSPS AND CORES 



for the Moore case. 

Fig. n shows an Aitoff projection of the angular posi- 
tion of these LG galaxies in Galactic coordinates. The 
size of each point is proportional to the 7-ray flux from 
Moore models, measured within a viewing angle of 1°. 
As shown by flSj, high energy 7-ray flux from nearby 
extragalactic halos with a Moore density profile can be 
larger than the 7-ray foreground produced by neutralino 
annihilation in our Galaxy. This is not an entirely triv- 
ial point as distant extragalactic objects are not resolved 
within the typical 1° angular resolution of the 7-rays de- 
tectors. To check whether our conclusions still hold for 
shallower density profiles, we have computed the 7-ray 
fiux using both Moore and NFW models for our LG halos 
and compared it to the Galactic foreground. In this work 
all fluxes are computed above 100 GeV within an angle 
An = IQ-^sr. Note that none of the particle physics 
assumptions affect the ratio between Galactic and extra- 
galactic fluxes that we will focus on. 



The previous results do not account for those astro- 
physical processes which may modify the density near 
the halo center, like the presence of a central SMBH. 
The effect of a SMBH on the central shape of a DM pro- 
flle is, in fact, rather controversial. The adiabatic growth 
of a SMBH of mass M, at the center of a DM halo would 
steep the slope of the density profile to ptaio , 
2.2 < 7 < 2.5 within a radius rh = GMt/a^, where a,, is 
the 1-D RMS velocity dispersion of the DM particles [33 . 
As a consequence, the expected 7-ray emission would be 
larger in galaxies hosting SMBHs. However, when one 
takes into account dynamical processes such as merger 
and orbital decay of an off-centered SMBH this 
central 'spike' disappears or is greatly weakened. More- 
over, in a CDM cosmology the assembly of black holes 
at the center of galaxies is the end-product of the hier- 
archical build up of their parent halos. In this scenario 
pregalactic black holes become incorporated through se- 
ries of mergers into larger and larger halos, sink to the 
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FIG. 2: Integrated 7-ray flux within AO, = 10~^sr as a func- 
tion of the angular distance from the GC. Filled triangles and 
flUed dots represent the fluxes from brightest extragalactic 
objects predicted assuming a Moore and a NFW proflle, re- 
spectively. The continuous and the dashed curves show the 
predicted Galactic foreground for the two model density pro- 
flies. 



center owing to dynamical friction, form a binary sys- 
tem Is^ and eventually coalesce ^si^ with the emission 
of gravitational waves. The process of formation and de- 
cay of black hole binaries transfers angular momentum to 
the DM particles, ejects them via gravitational slingshot 
effect and thus decreases the mass density around the bi- 
nary [4^ l4l| . The net result is either a constant density 
core '42I] or a shallower DM density profile ~ ^,-0.5 
within a radius rcut = 10 — 100 pc, with a consequent 
decrease of the neutralinos annihilation rate 14311. 



Observations indicate that SMBHs are only found in 
large galaxies. Indeed, among the LG galaxies only the 
Milky Way and M31 harbor a SMBH, while a possi- 
ble black hole in MSB would have a mass smaller than 
BOOOMq Q. Assuming that the presence of a SMBH 
creates a constant density core of similar size in both 
the dark and stellar components, we expect a small core 
Tcut = 0.38 pc in our Galaxy and one of rcut 10 pc 
at the center of M31 01 ■ a consequence, the SMBH 
phenomenon would only decrease the expected 7-ray flux 
from M31 by a factor of 3, still well above the Galactic 
foreground. In fact, our predictions do not change much 
even under the somewhat extreme hypothesis that all DM 
halos hosting LG galaxies have a constant density core 
of size 10 — 100 pc. Indeed, when repeating our calcu- 
lation using both NFW and Moore profiles truncated at 
10 and 100 pc we do find that the DM depletion near the 
halo center causes a significant decrease of the 7-ray flux 
which, however, does not prevent a few objects to shine 
prominently above the Galactic foreground. Results are 
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FIG. 3: Integrated 7-ray flux from the brightest LG members 
computed assuming a standard Moore proflle (filled triangles) 
and a Moore proflle truncated at 10 pc and 100 pc (flUed 
squares and filled dots, respectively). The continuous line 
show the expected Galactic annihilation foreground. 



shown in Fig. |3|in which the fluxes from Moore DM halos 
(flUed triangles), already shown in Fig. |21 are compared 
to the cases of Moore profiles truncated at rcut = 10 pc 
(filled squares) and at r^ut = 100 pc (filled circles). 



IV. HALO SUB-STRUCTURES 

Despite its remarkable success in explaining the large 
majority of current observations, the excess power on 
small scales in the ACDM concordance model has lead 
to what is sometimes called the CDM crisis. One aspect 
of this problem is the presence of a density spike at the 
center of DM halos, that is difficult to reconcile with the 
observed rotation curves of nearby dwarf galaxies [isj . 
A second aspect is related to the measured zero-point 
of the Tully-Fisher relation for spiral galaxies [46| . 
which is smaller than the CDM-predicted one. Whether 
these discrepancies are genuine or artifacts arising 
from observational biases (e.g. beam-smearing effect in 
determining the rotation velocities of dwarf galaxies) 
or numerical effects (e.g. limited resolution in CDM 
N-body experiments) is still matter of debate. A second 
problem, also related to the excess small-scale power in 
the CDM model, is the wealth of sub-galactic, virialized 
halos that do not appear to have an observational 
counterpart (the so called satellite catastrophe, [30|). 
Numerical simulations have indeed shown that a CDM 
halo of the size of our Galaxy should have about 50 
satellites with circular velocities Vdrr. > 20 kms^^ 
instead of the dozen observed ones |47|. Due to the 
hierarchical nature of the clustering process in a CDM 
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scenario, the mismatch is also found on larger systems 
of the size of the LG. Model predictions for halos with 
Vcirc > 35 kms~^ might be reconciled with observations 
by delaying the formation epoch of DM halos or by 
computing their peak rotation velocity with different 
techniques For lighter halos, however, the lack 

of a luminous counterparts could be either ascribed to 
astrophysical processes (supernovae explosion, stellar 
feedback, gas heating by ionizing radiation of cosmic 
origin) that suppress star formation or to a genuine 
disruption of the DM sub-halos by gravitational tides 
within their massive host. Finally, invoking Warm Dark 
Matter scenario decreases the small scale power and 
thus avoid the satellite catastrophe but would postpone 
the reionization of the universe to an epoch too late to 
be consistent with the recent WMAP observations of 
the cosmic microwave background |49l| . 

Here we account for the presence of a population of DM 
sub-halos and investigate their effect on the neutralino 
annihilation signal. Model predictions and related uncer- 
tainties are obtained using both numerical experiments 
and thoretical arguments. 

Large N-body simulations 

[13,113 

have shown that DM 
halos host a population of sub-halos characterized by a 
distribution function giving the probability of finding a 
sub- halo of mass m at distance r from the halo center, 
that can be modeled as follows: 1511: 



rish (r, m) = A 



-1.5 



(5) 



where r^'* is the core radius of the sub-halos distribution, 
niH is the mass of the parent halo and A is a constant set 
to obtain 500 sub-halos of mass > IO^Mq in a halo with 
TTiH = 2 • IO^^Mq [13 . The annihilation rate within sub- 
halos may contribute significantly to the total 7-ray fiux. 
According to 0, Q| the 7-ray Galactic foreground can 
be enhanced by over two orders of magnitude by the pres- 
ence of a population of sub-halos, an effect that may also 
help to explain the 7-ray flux excess found in EGRET 
data [13. To account for this effect, we have populated 
all the halos of our LG galaxies including our own with 
a population of sub-halos according to the distribution 
function |SJ| and computed their contribution to the 7- 
ray annihilation flux. For each sub-halo we use the same 
density profiles as its massive host. Current N-body sim- 
ulations do not allow to model Eq. lO very precisely. In 
particular, the mass clumped in sub-halos, rrici, and the 
value of r^'' are poorly constrained. To account for these 
model uncertainties we have generated several different 
sub-halos distributions that explore the edges of the pa- 
rameter space [mci,r^'^]. Following [13| we have set rrici 
equal to 10 % and 20 % of the total DM halo's mass, mn, 
and the core radius, rf^, to be 5 % and 30 % of the par- 
ent halo virial radius, r^^. Sub-halos were generated in 



a mass range [mmm, 0.1 mjy] and within a distance 
from the halo center. We set rrimin = 10^ Mq for Galactic 
halos and nimin = 1O^M0 for the the extragalactic ones 
since the slope of the sub-halo mass distribution, Eq. Q , 
is shallow enough to ensure that for a fixed value of md 
the contribution to the total 7-ray flux from halos with 
masses smaller than nrimin becomes negligible 0, [^ • 

Very high resolutions N-body experiments show that in 
a CDM scenario a large fraction of the sub-halo popula- 
tion survive the complete tidal disruption j23j . However, 
these sub-halos are tidally stripped of a fraction of their 
mass originating debris streams and possibly changing 
their original profile at all radii, thus lowerin g th e annihi- 
lation luminosity of the accreted systems [3SMI13|- be- 
cause of their finite resolution, the N-body experiments 
used to model Eq. ||SJ) underestimate the role of these 
tidal fields and therefore do not allow a reliable sampling 
of the low mass tail of the sub-halo distribution at small 
r. To better account for gravitational tides we follow 
[isj l and use the "tidal approximation", i.e. we assume 
that all the mass beyond the sub-halo tidal radius, rtid^ is 
lost in a single orbit without affecting its central density 
profile. The tidal radius is defined as the distance from 
the sub-halo center at which the tidal forces of the host 
potential equal the self-gravity of the sub-halo. In the 
Roche limit: 



rud{r) 



m 



(6) 



where r is the distance from the halo center. We trun- 
cate all sub-halos at their tidal radii and discard those 
with Ttid < for which the binding energy is positive 
and the system is dispersed by tides. Even in this case, 
we estimate the model uncertainties by considering three 
different cases: no tidal stripping, tidal radius computed 
at the mean orbital radius and tidal radius computed at 
the pericenter of the orbit (which is typically 1/5 of the 
mean orbital radius). 

The short-dashed histograms in the two bottom pan- 
els of Fig. ^ show the radial distribution of all Galactic 
sub-halos with M > 10^ Mq, distributed according to 
Eq. (O in which rud = 20% mn and r^'' = 5% for 
both a Moore (left) and a NFW (right) density profile. 
Including the tidal stripping effect destroys all sub-halos 
near the Galactic Center and modifies their radial distri- 
bution. The effect is larger for the Moore model (which 
is less concentrated than the NFW profile) and increases 
with the strength of the tidal field (i.e. when tides are 
computed at the pericenter of the orbit, long-dashed his- 
togram). The two top panels show the average tidal ra- 
dius, rtid of the sub- halo population that survive tidal 
stripping. 

The effect of sub-halos on the emitted annihilation flux 
is shown in Fig. 5a for the case of a Moore profile and in 
Fig. 5b for a NFW profile. The histogram drawn with 
a dotted line shows the diffuse 7-ray flux obtained by 
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FIG. 4: Number of Galactic halos (bottom) and their average 
tidal radius (top) as a function of the Galactic radius for a 
Moore (left) and NFW (right) density profile. Short-dashed 
histogram: no tidal stripping allowed. Continuous line: tidal 
radius computed at the mean orbital radius. Long dashed: 
tidal radius computed at the pericenter of the sub halo-orbit. 
In all cases shown nici — 20% mn and r"'^ — 0.05 r^^- 



adding the contribution to the smooth Galactic emission 
of a sub-halo population distributed according to Eq. Q 
in which rrid = 20% mn and rf^ = 5% r,^^. Filled dots 
show the flux from the most prominent external objects, 
computed by adding together the halo flux to that of all 
their sub-halos within 1° from the center and by neglect- 
ing gravitational tides. The presence of sub-halos dra- 
matically boosts up the level of the Galactic foreground 
(dot-dashed curve) by a factor 10-100, depending on the 
halo profile, while increasing the flux from extragalactic 
sources by a more modest factor 2-5. As a consequence, 
the Galactic foreground becomes dominated by the local 
sub-halo emission, resulting in the irregular 7-ray Galac- 
tic brightness profiles that outshine all extragalactic ob- 
jects but M31 and the LMC. 

These results, that are consistent with those of [l^ll^ . 
change significantly when tidal fields are accounted for. 
The net effect, that is more apparent for a Moore density 
profile, is that of decreasing significantly the Galactic 7- 
ray foreground, as shown by the histogram drawn with 
a continuous line in Figs. 5a and 5b. The flux from ex- 
tragalactic objects (filled squares) is comparatively less 
affected by tidal effects and in some cases may actually 
be enhanced when the chance of bright foreground extra- 
galactic sub-halos becomes non- negligible. As a result, 
taking into account the tidal stripping effects increases 
the prominence of extragalactic sources over the Galactic 
foreground. If one further accounts for the eccentricity 
of the sub-halo orbits and assume tidal stripping occur 
at their pericenter then the Galactic 7-ray foreground 
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FIG. 5: Integrated 7-ray flux within Afl = 10~^sr for Moore 
(top) and NFW (bottom) DM halos populated with sub-halos. 
Galactic (histograms) and extragalactic (symbols) fluxes are 
shown for the different sub-halo populations described in the 
text. 



(dashed histogram) and the fiux from extragalactic ob- 
jects (filled triangles) further decrease to a level similar 
to that of Fig. |21 

To better quantify the uncertainties in modeling the 
sub-halo population we list in Table ^ the fractional 
variation of the 7-ray flux from the extragalactic sources 
when varying the model parameters. The quantity 
displayed is 5$ = {^"'^ - <i>°)/$°, where is the flux 
from a smooth profile while the fiux also accounts 
for the sub-halo population. We have focused on the 
two brightest LG objects, M31 and the LMC. It is worth 
noticing that, due to their different distances from our 
Galaxy, while the fiux from M31 is produced in the 
central 16-36 Kpc (depending on the model profile), the 
LMC flux is produced in a much smaller region of radius 
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1-1.4 Kpc. The various models explored are listed in 
Table m and are characterized by a label in which the 
first two digits indicate the extent of the core radius of 
the sub-halo distribution (e.g. 30 means r*'* — 30% r^^), 
the subsequent two digits characterize the mass clumped 
in sub-halos (e.g. 20 for m^i = 20% m_f/) and the last 
two letters indicates the model tidal field (dO — no tidal 
effects, dd — tidal radius computed at the mean orbital 
radius, dh = tidal radius computed at the pericenter). 





5$ 


Galaxy 




0520dd 


iiOlOdd 


:i020dd 


0520cl0 


Orj20dr3 


LMC (Moore) 


-0.12 


-0.21 


-0.12 


-0.21 


+1.35 


-0.22 


LMC (NFW) 


-0.07 


-0.18 


-0.09 


-0.18 


+0.90 


-0.18 


M31 (Moore) 


-H0.25 


+0.52 


+0.03 


+0.03 


+2.69 


-0.10 


M31 (NFW) 


-hO.78 


+1.59 


+0.02 


+0.05 


+1.93 


+2.09 



TABLE II: The effect of including sub-halos on the 7-ray 
flux from LMC and M31 in different halo models. Each en- 
try specifies the fractional fiux variation with respect to the 
smooth halo model. The choice of parameters is encoded in 
each model label as specified in the text. 



experiments designed to detect the products of generic 
atmospheric showers. These detectors were not originally 
meant to explore exotic physics, yet recent advances in 
technology have allowed to reach unprecedented high 
sensitivities and they can now be used to detect high 
energy photons from neutralino annihilation. 

Cerenkov detectors. Electrons generated in showers 
emit Cerenkov light during their passage through the at- 
mosphere. It reaches the ground in the form of a ~ 100 
meters wide light pool. Cerenkov detectors have good an- 
gular and energetic resolution, operate with a duty cycle 
of about 10% and have a very good hadronic background 
rejection. Their small angle of view (~ 5°) makes them 
suitable for the observations of point sources like LG 
galaxies. Currently operating Cerenkov telescopes such 
WHIPPLE H] and HEGRA do not have enough 
sensitivity to detect faint extragalactic sources. Planned 
Active Cerenkov Telescopes (ACT) such as HESS [s^, 
MAGIC m and VERITAS ^ with their lower en- 
ergy threshold (~ 50 GeV) and greater sensitivities make 
them suitable for our purposes. 



In all dd models considered, the flux from M31 is in- 
creased by the presence of sub-halos while the flux from 
LMC decreases. This discrepancy reflects the different 
effect of the tidal field that, as we have seen, disrupts 
all the sub-halos near the halo center and thus decreases 
the flux from the nearby LMC, while the flux from M31, 
which is produced in a broader region, receives a signifi- 
cant contribution from the nearby, foreground sub-halos 
in M31. Varying rrici in the allowed range does not change 
the fluxes appreciably, while increasing the core radius of 
the sub-halo distribution decreases the M31 flux signifi- 
cantly but does not change significantly the flux received 
from the nearer LMC. Finally, neglecting the tidal field 
(do model) has the rather obvious effect of increasing the 
flux from both extragalactic objects. Enhancing the tidal 
disruption effect (d5 model) generally decrease the anni- 
hilation flux, apart from the case of M31-NFW, in which 
the flux is enhanced by the combined effect of the halo 
distance and the large DM concentration. 

V. DETECTABILITY 

In this section we explore the actual detectability of 
the 7-ray fluxes from both the extragalactic sources and 
the Galactic foreground. 

Satellite borne experiments, such as GLAST JSJ] , have 
very small effective areas at energies > 250 GeV and 
thus are inefficient in revealing high energy photons 
from neutralino annihilation. The need for larger 
collecting areas can only be met by those ground-based 



Large Geld of view arrays. Unconventional air shower 
arrays at energies < TeV constitute viable alternative to 
Cerenkov telescopes. Unlike standard air shower arrays, 
these detectors can reconstruct lower energy showers 
by dense sampling of the shower particles. This can be 
achieved by operating the detector at very high altitude 
in order to approach the maximum size development of 
low energy showers and by using a full coverage layer of 
counters. Operating large field of view arrays like ARGO 
[Hfij i and MILAGRO have worse angular resolution 
and background rejection than ACT detectors, but are 
able to explore the whole sky with a duty cycle of about 
100 %. 



Due to the complementarity of ACTs and large field of 
view arrays we explore and compare the capabilities of 
these detectors of revealing the 7-ray flux from neutralino 
annihilation in extragalactic objects [s^l- AH analyses, 
including our own, have shown that the best place to 
look for DM annihilation signals is the GC, and indeed 
presently available and future detectors should be able to 
detect such 7-ray emission. However, most of them are 
(or will be) located in the Northern hemisphere where 
the GC can only be seen for short periods and at uncom- 
fortably large zenith angles. In Fig. Elthe angular posi- 
tions of our extragalactic objects (filled dots) are shown 
in Galactic coordinates, superimposed to the areas that 
can be seen with a zenith angle smaller than 30° from var- 
ious experimental sites (shaded areas). In the following 
we estimate the detectability of 7-ray fluxes from either 
external LG galaxies or sub-halos within our Galaxy. 
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FIG. 6: Angular positions of our 44 LG galaxies (filled dots) 
inclusive of M87 and the GC in galactic coordinates. Shaded 
regions have zenith angles smaller than 30° at the ARGO 
site (upper panel), HESS site (lower left) and an hypothetical 
Antarctic site (lower right). Terrestrial longitude and latitude 
of each sites are indicated above each panel. 



Detectability of external galaxies 

In the following we will focus on the detectability of 
M31 since it is the extragalactic object with the largest 
7-ray flux visible from the Northern hemisphere. 

The detectability of the continuum flux from M31 is 
computed by comparing the number of 7 events expected 
from this source to the fluctuations of background events. 
This ratio, which we call a in the following, is given by: 



dE 



-dEde 



Jbck ^bcki^' 



(7) 



.A 



dE 



•--dEde 



where Ts is the time during which the source is seen 
with zenith angle 9 < 30° and eAsi — 0.7 is the fraction 
of signal events within the optimal solid angle Afl cor- 
responding to the angular resolution of the instrument. 
The effective detection areas A'^-^^ for electromagnetic 
and hadronic induced showers are defined as the detec- 
tion efflciency times the geometrical detection area. For 
the case of a large field of view array we have considered 
A^-f-l" = 10^ — 10^ m?, depending on the 7-photon 
energy, while for a Cerenkov apparatus we have assumed 
a conservative effective area of 10'' rri^. Note that 
while the latter can be increased by adding together 
more Cerenkov telescopes, the former is intrinsically 
limited by the number of hits reaching the ground and 
cannot be much greater than our fiducial value. Finally 
we have assumed an angular resolution ((> = 1° and 
hadron-photon identification efficiencies of e = 75% for 
the case of a large field of view array and 4> = 0.1° and 



e = 99% for the Cerenkov apparatus. These values are 
appropriate for the energy range of interest. 

Backgrounds. In our analysis the M31 flux is much 
smaller than electron, hadron and diffuse 7-ray back- 
grounds that must then be taken into account. We have 
considered the following values for the background levels: 



d(t>^ 



IA9E~'^'^^- 



P 



dQdE cm^ s sr 

for the hadronic background [H^ . 

J J. p. 

6.9 X IQ-^E-^-^ ^ 



GeV 



dQ.dE 

for the electron background 'b^I 



CTO^ s sr 



GeV 



and 



'^dif fuse 

dfldE 



1.38 X 10""^; 



-2.1 
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crn^ s sr GeV 



(8) 



(9) 



(10) 



for the diffuse extragalactic gamma emission, as extrap- 
olated from EGRET data at lower energies [65j. 
The Galactic diffuse emission along M31 line of sight 
[li j turned out to be negligible with respect to the 
annihilation signal and has been neglected. 
In Fig. [3 the various backgrounds (dashed lines) are 
compared to the integrated energy spectra of 7-photons 
from DM annihilation in M31, computed assuming 
a Moore DM profile and a population of sub-halos 
unaffected by tidal effects. The various curves refer to 
different neutralino masses, as specified by the labels. 
Filled dots represent the expected 7-line fluxes. 

Fig. IHl compares the expected integrated 7 flux from 
M31 (dot-dashed curve) to the 5 a sensitivity curves for 
the large field of view array (upper line) and the Cerenkov 
telescope (lower line). The 5 a detection curves assume 
1 year of data taking for the large field of view array and 
20 days pointing for the Cerenkov telescope. 
The results show that annihilation signals from M31, 
which are too faint for large field of view arrays, are well 
within reach of next generation Cerenkov telescopes. 



Detectability of sub-halos in our Galaxy 

The only way of detecting signals from annihilation in 
Galactic sub-halos is through blind searches since model 
predictions specify the sub-halo distribution but not their 
precise spatial positions. With this respect, the use of a 
large field of view array has the advantage of collecting 
signals from a wide area of the sky, while ACTs would 
only allow to survey a much smaller region. 

To estimate the chances of detecting a 7-ray signal 
with cither detectors we performed ideal observations 
with both instruments and computed the significance of 
an annihilation signal using Eq. Q. 



9 




60 70 80 90100 200 300 400 500 600 700 

Energy (GeV) 



FIG. 7: Integrated spectra of photons from DM annihilation 
in M31 for different values of the neutralino mass, specified by 
the labels, computed assuming a a Moore profile filled with 
a population of sub-galactic halos unaffected by tides. Both 
the continuum emission (curves) and the 7-line contribution 
(filled dots) are shown. The relevant backgrounds (dashed 
lines) are shown for comparison. 
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FIG. 8: 5a sensitivity curves for a high altitude full coverage 
air shower array (upper line) and for a Cerenkov cell (lower 
line), computed using the parameters specified in the text. 
The curve shows the same fiux from M31 displayed in Fig. 6 
for the case of a neutralino mass of 1 Tev. 

We restrict our observations to all directions with 
\b\ < 30° if 1^1 < 40° around the GC and |6| < 10° 
elsewhere provided that they are within a zenith angle 
of 30°. These constraints allow us to ignore the Galactic 
diffuse photon contribution to the background emission 
and the effects of the atmospheric depth. 

For the case of a large field of view array we have 



chosen the ARGO experimental site (6.06 h E, 30.18° 
N) in Tibet and considered the case of a 365 days 
observation with a detector having an angular resolution 
of 1°. Each point in Fig. |51 represents the significance 
of a 7-ray signal detection in such experiment, plotted 
as a function of ip, the angle from the GC. Each filled 
dot accounts for the annihilation flux received from all 
Galactic sub-halos with the same ip that have entered 
the detector field of view during the observation. The 
Galactic signal in absence of sub-halos is also plotted for 
reference (open dots). The significance detection level is 
very small at all ip but in a couple of cases where a very 
nearby sub-halo crosses the line of sight. The results 
show that it would be impossible, using a large field of 
view array, to identify the Galactic sub-halo origin of 
any detected 7-ray flux. 

For the case of an ACT we have chosen the VERITAS 
site (8.87 h W, 33° N) in Arizona. In this case we have 
considered the results of a 20 days "stacking" observation 
in which an instrument with an angular resolution of 0.1 
degrees points toward different directions characterized 
by the same value of "0- The results of this experiment 
are represented by the triangles in Fig. |51 Since each 
point represents a 20 days observations we have sampled 
only a few angles ip. Although the significance of the 
ACT detections is higher than in the case of a large 
field of view array, the serendipitous nature of the 7-ray 
sources combined with the small extension of the sur- 
veyed area make it very hard to detect the annihilation 
signal from Galactic sub-halos using Cerenkov telescopes. 

We then conclude that only satellite borne experiments 
like GLAST will have the chance to detect 7-ray photons 
from Galactic sub-halos and to study their spatial distri- 
bution 

VI. DISCUSSION AND CONCLUSIONS 

In this work we have computed the expected 7-ray flux 
from DM annihilation in the nearest LG galaxies, under 
the assumption that the bulk of DM is constituted by 
neutralinos. Predicted fluxes depend on a number of as- 
sumptions like the shape of the DM halo density proflles, 
the possible existence and the distribution of a popula- 
tion of sub-galactic halos and the dynamical influence 
of SMBHs located at the halo centers. Uncertainties in 
these theoretical models that propagates into flux pre- 
dictions have been evaluated by systematically exploring 
the space of the currently accepted model parameters. 

It turns out that, in all cases considered, the fluxes 
from the Small and the Large Magellanic Clouds, M31 
and M87 are well above the level of the 7-ray annihila- 
tion Galactic foreground. Uncertainties in the predicted 
Galactic emission are quite large (they vary by a fac- 
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FIG. 9; Significance of the 7-ray signal from Galactic sub- 
halos as a function of the angle of view from the GC, for a 
large field of view array (dots) and for a Cerenkov telescope 
(triangles). The considered sub-halo distribution is not af- 
fected by any tidal effect. 



stellar dynamics and strong lensing data in a sample 
of brightest cluster galaxies jjT^ , and from microlensing 
events toward the GC (a = 0.4, in the hypothesis that 
the MW halo is spherically symmetric |7J|). However, 
estimates based on weak gravitational lensing in X-ray 
luminous clusters indicate a much steeper profile with 
0.9 < a < 1.6 [73. 

The previous examples show that the current observa- 
tions neither support nor reject profiles as shallow as 
a = 0.5 or 0. Therefore, while we regard the NFW 
and the Moore profiles as the most likely cases based 
on N-body results, we also want to keep an eye on the 
shallow profile case. For this purpose we have estimated 
the expected 7-ray flux from both M31 and LMC in the 
case of shallow profiles. If a = 0.5 then the M31 flux is 
r^ 50 time smaller than in the NFW case and decreases 
by a further factor of 10 with a flat core a = 0. Results 
for the LMC are very similar. Therefore, if the slope 
of the density profile would turn out to be as shallow 
as a = 0.5 or more, then no ground-based experiments 
would have the chance of detecting extragalactic 7-ray 
annihilation fluxes. 



tor ^ 100) and mainly result from the poorly known 
modeling of the DM sub-structures within our Galaxy. 
This scatter actually increases when considering more ex- 
otic possibilities that would further boost up the Galac- 
tic foreground like the existence of caustics in the DM 
distribution or the presence of a population of mini 
sub-halos allowed by the very small cutoff mass in the 
CDM spectra Neutralino annihilation in sub-halos 
contributes less to 7-ray flux from extragalactic sources, 
resulting in a much smaller (by a factor of ^ 10) model 
uncertainties and more robust flux predictions. 

In section II we have assumed that the inner slope 
of the DM density profile p{r) cx r"" is in the range 
1.0 < a < 1.5 following the indications of numerical 
experiments rather than using observational constraints. 
The reason was that current observations do not provide 
a self-consistent scenario. Indeed, recent determina- 
tions based on the X-ray properties of the intercluster 
medium seem to indicate a steep DM density profile 
1.0 < a < 2.0, provided that clusters with disturbed 
X-ra y su rface brightness are removed from the sample 
[Hil l69l |. On the other hand, the rotation curves 
measured in a sample of LSBs by 0| are consistent, 
on average, with a shallower density profile a = 0.2. 
However, when the full extent of the rotation curves is 
taken into account rather than the inner region only, 
then the recovered densit y p rofiles lie between the 
NFW and the Moore ones |71|. Another study of high 
resolution Ha rotation curves for dwarf and LSB |23| 
has ruled out neither a profiles as steep as a = 1 nor a 
constant density core. Finally, indications for shallow 
profiles {a = 0.52 ± 0.3) were found by combining 



In the hypothesis of both a NFW and a Moore DM 
density profile we have estimated the possibility of de- 
tecting the annihilation signals from Galactic sub-halos 
and extragalactic sources using those ground based in- 
struments sensitive enough to detect 7-photons in the 
^ 50 GeV - few TeV energy band, such as ACTs or 
large field of view arrays. We have found that the ex- 
pected fluxes from Galactic sub-halos are too faint to 
be detected by ground-based observatories. Their exis- 
tence and distribution could instead be probed by next 
generation satellite borne experiments such as GLAST. 
On the other hand, ground-based experiments should be 
able to detect a few extragalactic objects such as LMC 
and M31. In particular, we have been focusing on M31 
which turned out to be the brightest extragalactic object 
visible in most experimental sites. We have shown that 
next generation Cerenkov telescopes should be capable of 
revealing the 7-ray flux from M31 at a significance level 
~ 5(7 in a 20 days pointing observation. These results 
have been obtained using a quite conservative observa- 
tional set up and theoretical modeling. Increasing the 
number of Cerenkov cells (as already planned for some 
future experiments) or the serendipitous superposition 
of a nearby Galactic sub-halo would certainly increase 
the chance of detecting an annihilation signal along the 
M31 direction. The capability of detecting extragalactic 
signals is not shared by large shower arrays, which due 
to their lower sensitivity will not be able to detect such 
a faint source. It is worth stressing that the GC does 
remain the best place to consider for detecting DM an- 
nihilation signatures. Unfortunately, and unlike M31, it 
is not visible from the Northern hemisphere where most 
of the ground-based detectors are (or will be) located. 
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Finally, it is worth stressing that our predictions as- 
sume that DM is made of neutralinos. However, the 
natural factorization of particle physics and cosmology 
in the problem addressed in this paper makes straight- 
forward to extend our approach to other kinds of weakly 
interacting DM candidates. 
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